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Q\ Abstract 
O 

In order to describe correctly the interplay of extrinsic and intrinsic spin-orbit mechanisms to the spin Hall effect, it is necessary 
to consider different sources of spin relaxation. We take into account the spin relaxation time top due to the Dyakonov-Perel 

Q^mechanism as well as the Elliot- Yafet spin-relaxation time t s due to the spin-orbit scattering from impurities. The total spin Hall 

(D conductivity depends crucially on the ratio t s /tdp- 
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1. Introduction 

In a recent paper[ 1 1, we have provided a general framework 
to describe the electric-field control of the electron spin in a 
two-dimensional electron gas (2DEG) with diffusive electron 
spin dynamics when both intrinsic and extrinsic spin-orbit inter- 
actions are present. The extrinsic mechanism for the spin Hall 
effect arises to first order in the spin-orbit coupling strength, /& 
(see below)||2][3|. However, when also the intrinsic mechanism 
is present|4] [5] |6), a first-order calculation is no longer suffi- 
cient, since it is important to consider the Elliot- Yafet spin re- 
laxation due to the spin-orbit scattering from impurities, which 
arises in the order A A Q . Here, by using a Keldysh Green function 
approach[7|, we provide a microscopic basis of the equation 
of motion for the spin density that governs the interplay of in- 
trinsic and extrinsic mechanisms. The layout of the paper is as 
follows. In the next section we introduce the model Hamilto- 
nian and the method. In section 3 we derive an expression for 
the spin current and an associated continuity equation. In sec- 
tion 4, we consider the specific case of the spin Hall effect and 
derive a formula for the spin Hall conductivity. Finally, we state 
our conclusions in section 5. 



• 2. The model and the method 

X 

^ In the presence of both extrinsic and intrinsic spin-orbit in- 
teraction as well as normal potential scattering from impurities 
V(x), the Hamiltonian for the 2DEG can be written in terms of 
a spin-dependent [S U(2)] vector potential A 



H 



2m 



A p 



V(x). 



(1) 



The vector potential is the sum of intrinsic (Rashba type) and 
extrinsic contributions 



A = m cr x 



tre 



A 2 

+ ~jd x V(x) 



= -AV, 
2 



(2) 



where cr is the vector of Pauli matrices and a = x,y,z. The 
coupling constants a and A% characterize the strength of the 
intrinsic and extrinsic spin-orbit interaction, respectively. The 
advantage of introducing the S U(2) vector potential is that one 
can immediately derive a continuity equation for the spin den- 
sity 

8,s a + d x -T +e abc k b -f = 0, (3) 



where the spin density is defined in terms of the Heisenberg 
field operators s a = (l/2)(i/r(x, f)cr°i/r(x, f)) and the spin current 
has the expression j" = (l/4)(i^'(x, f){o _fl ,v}(/'(x, f)), with the 
velocity operator given by 



p-A 



(4) 



Although Eq. ([3]l is formally exact, it cannot directly convey in- 
formation for the disorder averaged spin current, since the dis- 
order potential appears also in the vector potential A. In order 
to carry out the average over the disorder, it is convenient to use 
the Green function approach as developed in Ref. Q. It is also 
useful to explicitly separate the spin-orbit interaction due to the 
intrinsic mechanism from that due to the extrinsic mechanism. 
To this end we define a space-independent vector potential 



A = am cr x e, 
and a spin-dependent disorder potential 



(5) 
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U(x) = V(x) -jcrx d x V(x) ■ p 



(6) 
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Figure 1: Lowest order diagrams for the disorder averaged self-energy, (a) 
Born approximation, (b) Third order contribution. The full dot indicates the 
insertion of the potential U(x) (cf. Eq.JSJ), whereas the dashed line represents 
the disorder average. 



This contribution leads to the standard elastic scattering time 
r -1 = 2nNonivl |9|, where No is the single-particle density of 



states of the Fermi gas. To first order in Ai, Eq. (Ill yields the 
three terms 

A 2 

^a(P.X) = -i^-nivle abc ^[p a o- h p' c ,G(p',x)], (13) 



^(P,x) 



A. 

= -fn i v 2 e abc Y u {o-bP' c ,d a G($',x)}, (14) 



The disorder averaged spin density is given by 
2 J 2tt J (2n) 2 



s a (x, t) 



1 f de C d 2 ] 

2 J 2^J (2* 



Tr(cr fl G < (e,p,x,f)), (7) 



where G < (e, p, x, f) is the lesser component of the disorder aver- 
aged Green function, here given in Wigner coordinates! 8]. The 
current density is written as a sum of two terms, j a = jJJ + j° r , 
associated to the normal, vo = (p - A)/m, and to the anomalous 
disorder-dependent components of the velocity, respectively. 
Explicitly the normal current is 



j ° (X '° = 4 /£ /^TrK { v„,G< (e ,p,x, f )}], 



(8) 



where the symbol {, } denotes the anticommutator. Since the 
anomalous velocity contribution to the current contains explic- 
itly the disorder potential, its expression can be obtained only 
after specifying the disorder model and the approximations 
used. Later on, we will identify the expression for the anoma- 
lous current from the self-energy terms contributing to the con- 
tinuity equation. 

In Wigner coordinates and after a gradient expansion, the 
equation of motion for the Green function G(e, p, x, t) reads 



i<9,G + id x ■ { 



P-A 

2m 



G + 



A p 



[%G]. 



(9) 



The Green function has a matrix structure in both the Keldysh 
and spin space, and the symbol [, ] indicates a commutator. An 
external electric field E is included via the substitution d x — > 
<9 X - eE<9 £ . The self-energy S entering the RHS of the Eq.(Mi 
contains all the effect of the disorder potential defined in Eq.d6]i. 
We assume the standard model of uncorrelated impurities with 



V(x 1 )V(x 2 ) = n i v^<5(x 1 -x 2 ), 



(10) 



where is the impurity concentration and Vo the scattering am- 
plitude. In the Born approximation it is sufficient to define the 
moments of the disorder distribution up to the second order as 
done in Eq.([T0|). At the level of the Born approximation, the 
self-energy is given by the diagram of Fig. [TJa) and reads 



S(xi,x 2 ) = U(x l )G(x 1 ,x 2 )U(x 2 ). 
To zero order in A 2 y Eq. (Hi yields 



i°(p,x)=n i v 2 Y J G(p',x). 



(11) 



(12) 



A 2 

£i(p,x) = -jn iV 2 e abc Y}PaO-b,d c G{ 9 ' ,x))- (15) 
p' 

again the external electric field is included in the space deriva- 
tive d a = d Xa -eE a d t . The self-energy t} a is related to the mech- 
anism for "swapping of spin currents" HOI . Si and E' are side- 
jump contributions. Notice that to lowest order in the gradi- 
ent expansion, only the "swapping term" is different from zero. 
The side-jump contributions arise when considering the next- 
to-leading order in the gradient expansion. Finally, to second 
order in A^, Eq. (Hi yields 



t 2 



16 



^ JVG(k,x)<r%xk)2, 

1 k 



(16) 



which corresponds to the Elliot- Yafet spin relaxation mecha- 
nism. 

Skew-scattering arises beyond the Born approximation J3] [3), 
starting from the third order diagram shown in Fig. |TJb). By 
using a disorder model with third moments different from zero 
and defined by[3 | 



V(xi)V(x 2 )V(x 3 ) = «iv^(xi - x 2 )<5(x 2 - x 3 ), (17) 
the diagram of Fig.jTJb) yields to first order in A 2 r the terms 

^2 

^ss,a = i>^G(k,x)G(k',x)pxk'-(r (18) 

k,k' 



f 1 

^SS.b 



J/ Ii v^kxp-o-G(k,x)G(k',x). (19) 



k.k' 



To make contact with the diagrammatic language of the Kubo 
formula, we notice that the two self-energies ( T8~fl9 i corre- 
spond to the diagrams of Fig. 2 in Ref.j3). 

3. The continuity equation and the spin current 



The insertion of Eqs. ( fT2|fT6[ i and ( |18]{T9[ l into the equation- 
of -motion (|9]l allows to derive a continuity equation for the spin 
density polarized along the a-axis. After integrating over (e, p) 
and taking the trace of Eq. d9), one obtains 



i 



(20) 



The gradient together with the commutator on the LHS of 
Eq. (|9]l are the origin of the covariant derivative of the spin cur- 
rent 



(21) 



Most of the contributions to the self-energy disappear in the in- 
tegrated equation. The only terms surviving the integration are 
those originating from and 2 2 . The term containing l, l b of 
Eq. ( 14 1 yields a contribution which can be written as a diver- 
gence and hence, apart from a minus sign, defines the anoma- 
lous contribution to the current j" ]v . The spin relaxation time, 



which is obtained from Eq. (16 1, only applies to the in-plane 
spin components (cf. the two cr z matrices before and after the 
Green function) and reads r^ 1 = T~ l (Aopf/2) 4 . 

Clearly, also the explicit expression for the spin current can 
be obtained starting from Eq. |9]). To do this, we followed the 
procedure of Ref. [7|, namely we integrated Eq. |9]) over £ = 
p 2 /2m - p in order to obtain a Boltzmann-like equation for the 
quasiclassical Green function, 

|(p,x) = l ~ J d^G R/A (p, 



x). 



(22) 



Technically the equations simplify considerably by using that 
the retarded and advanced quasiclassical Green functions are 
given by a constant, g R,A m ± 1 . . 

After some more steps we found, in the diffusive limit, when 
the spin splitting due to the intrinsic spin-orbit interaction, ap F , 
is smaller than the disorder broadening, r _1 the spin current as 



Ji 



Dd iS " + cr£ t s iab E b + y a u Ej. 



The contribution from the extrinsic spin-orbit mechanism has 
the form first predicted by Dyakonov and Perel[ll | and is due 
to crf^ = cr^f + cr s s H with the standard expressions for the side- 
jump (sj) and skew-scattering (ss) contributions 



4 er' 



1 



-(p F l)(27TN v )o- 



sH 



(24) 



cr D = 2e 2 N(,D being the Drude conductivity and D = v 2 f t/2 
the diffusion coefficient. It is worthwhile pointing out that the 
anomalous disorder-dependent contribution to the current is one 
half of the side-jump contribution, f avi = (l/2)cr^e, fli £' i . 

The intrinsic spin-orbit coupling yields a contribution to the 
spin current due to the SU(2) magnetic field B fl = d x x A", 



Tcr D 
8 me 



(B fl x E);. 



(25) 
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Figure 2: Spin Hall conductivity as a function of the strength of the in- 
trinsic spin-orbit coupling. The numbers are obtained for GaAs with den- 
sity n = 10 12 /cm 2 , mobility fi = 10 3 cm 2 Vs, extrinsic spin-orbit coupling 
Aq = 4.7 X 10~ 8 cm and assuming scattering from positively charged impurities 
with Nqvq = -1/2. The spin Hall conductivity in the absence of the Rashba 
coupling is cr sH (0) = -3 X l(r 7 /n. 



where the first term originates from the covariant derivative de- 
The presence of the in-plane spin density, 



fined in Eq.(21 



requires that Eq. ( 26 1 must be solved together with the y- 



(23) component of the continuity equation ( 20 1 



d t s y + 2ma\jy 



1 

+ — s 

T 



0. 



(27) 



In the absence of the extrinsic mechanism (Ao = 0), as it has 
been noticed previously fl6l[T7l[T8 l. the above equation implies 
the vanishing of the spin current and spin Hall conductivity! 19, 
|20l |2T| [22) . As a result the in-plane spin polarization acquires 
the electric-field dependent value predicted by Edelstein|23 1. 
In the presence of both the extrinsic and intrinsic mecha- 



nisms, by solving Eqs. (26 1 and (27 1 together, one gets j~ 
a sH E with 



a sH = 



l/r s 



1/t s + 1/tdp 



<+^?4<K<' (28) 



where we have introduced the Dyakonov-Perel relaxation time 



For the Rashba model the only non-zero component of B" is t dp = (2ma) 2 D. To first order in A 2 , r; 1 = and one obtains 



B\ = 2(2ma) 2 leading to yfjEj = s iab o-^E b with cr^ = 
j^{apff) 2 . The extra terms of Ref. 11101 . associated to t} a , are 
not relevant in the present context and therefore have been ig- 
nored in Eq. (j23j 



4. The spin Hall conductivity 



After these formal aspects on how to derive Eqs. (20i and 



( 23 1 we now apply the equations to a specific problem, namely 
the spin Hall effect 0J] [Til US d [Bl . More precisely we study 
the transverse spin current that is generated by a uniform elec- 



tric field applied along the x-axis. According to Eq.(23 i the spin 
current flowing along the y-axis reads 



2maDs y + ((j% t + a£ t )E, 



(26) 



<t sH = (l/2)cr s s H for any a + 0, in agreement with |4]. Clearly, 
the ratio tJtdp acts as a control knob of the spin Hall con- 
ductivity. This suggests that by adjusting the constant a by a 
suitably applied gate voltage, one can vary experimentally the 
magnitude of the spin Hall current. For illustrative purposes in 
Fig. |2j we plot the spin Hall conductivity for GaAs as a function 
of the dimensionless parameter 2ap^T. 

5. Conclusions 

In summary we have presented a microscopic derivation of 
the equation of motion for the Keldysh Green function in the 
presence of both extrinsic and intrinsic spin-orbit interaction as 
well as scattering from impurities. In particular we have given 
explicit expressions for the disorder averaged self energy at the 



3 



level of Born approximation and to the first order beyond the 
Born approximation. This has allowed us to derive an expres- 
sion for the spin current and the associated continuity equa- 
tion. It has been shown that the ratio of the Elliot- Yafet and 
Dyakonov-Perel spin relaxation times is the important param- 
eter controlling the interplay of extrinsic and intrinsic mecha- 
nisms to the spin Hall effect in a 2DEG. 
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